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Abstract 

Recently Vladimir Novikov found a new integrable analogue of the 
Camassa-Holm equation, admitting peaked soliton {peakon) solutions, 
which has nonlinear terms that are cubic, rather than quadratic. In this 
paper, the explicit formulas for multipeakon solutions of Novikov's cubi- 
cally nonlinear equation are calculated, using the matrix Lax pair found 
by Hone and Wang. By a transformation of Liouville type, the associated 
spectral problem is related to a cubic string equation, which is dual to 
the cubic string that was previously found in the work of Lundmark and 
Szmigielski on the multipeakons of the Degasperis-Procesi equation. 



1 Introduction 

Integrable PDEs with nonsmooth solutions have attracted much attention in 
recent years, since the discovery of the Camassa-Holm shallow water wave 
equation and its peak-shaped soliton solutions called peakons [5]. Our pur- 
pose in this paper is to expHcitly compute the multipeakon solutions of a new 
integrable PDE, equation l|3.ip below, which is of the Camassa-Holm form 
Ut — Uxxt = F{u,Ux,Uxx, ■ ■ ■), but has cubically nonlinear terms instead of 
quadratic. This equation was found by Vladimir Novikov, and published in a 
recent paper by Hone and Wang jl9i |. 

We will apply inverse spectral methods. The spatial equation in the Lax 
pair for Novikov's equation turns out to be equivalent to what we call the dual 
cubic string^ a spectral problem closely related to the cubic string that was 
used for finding the multipeakon solutions to the Degasperis-Procesi equation 
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\27 \ [28 l [23]. Once this relation is established, the Novikov peakon solution can 
be derived in a straightforward way using the results obtained in [23]. The 
constants of motion have a more complicated structure than in the Camassa- 
Holm and Degasperis-Procesi cases, and the study of this gives as an interesting 
by-product a combinatorial identity concerning the sum of all minors in a sym- 
metric matrix, which we have dubbed the Canada Day Theorem (Theorem l4.lt 
proved in Appendix [Aj . 

The peakon problem for Novikov's equation presents in addition one impor- 
tant challenge. Unlike its Camassa-Holm or Degasperis-Procesi counterparts, 
the Lax pair for the Novikov equation is originally ill-defined in the peakon sec- 
tor. The problem is caused by terms which involve multiplication of a singular 
measure by a discontinuous function. We prove in Appendix [B] that there ex- 
ists a regularization of the Lax pair which preserves integrability of the peakon 
sector, thus allowing us to use spectral and inverse spectral methods to obtain 
the multipeakon solutions to the Novikov equation. This regularization problem 
has a subtle but nevertheless real impact on the formulas. In general, the use of 
Lax pairs to construct distributional solutions to nonlinear equations which are 
Lax integrable in the smooth sector but may not be so in the whole non-smooth 
sector is relatively uncharted territory, and the case of Novikov's equation may 
provide some relevant insight in this regard. 

2 Background 

The main example of a PDE admitting peaked solitons is the family 

Ut - Uxxt + {b+ l)uUx = bUxUxx + UUxxx, (2.1) 

often written as 

mt + nixU + bmux = 0, m = u — Uxx, (2-2) 

which was introduced by Degasperis, Holm and Hone [10], and is Hamiltonian 
for all values of b ^7\. It includes the Camassa-Holm equation as the case 6 = 2, 
and another integrable PDE called the Degasperis-Procesi equation [TTl [lO] as 
the case 6 = 3. These are the only values of 6 for which the equation is integrable, 
according to a variety of integrability tests [TT | ISO t fTS t 120]. (However, we note 
that the case 6 = is excluded from the aforementioned integrability tests; 
yet this case provides a regularization of the inviscid Burgers equation that is 
Hamiltonian and has classical solutions globally in time [4].) Multipeakons are 
weak solutions of the form 

n 

u(a;,i) =51™'^^"''^"''''*^'' (2.3) 

i=l 

formed through superposition of n peakons (peaked solitons of the shape e^'^^l). 
This ansatz satisfies the PDE (|2.2p if and only if the positions (a;i, . . . , x„) and 
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momenta (mi, . . . , m^) of the peakons obey the following system of 2n ODEs: 

n n 

ifc = ^m,e-l'"'-"^*l, Thk =^ {b-l)mk^m, sgn{xk~x^)e-^'"'-''^^. (2.4) 
1=1 1=1 

Here, sgnx denotes the signum function, which is +1, —1 or depending on 
whether x is positive, negative or zero. In shorthand notation, with {f{x)) 
denoting the average of the left and right hmits, 

(/(^)) = ^(/(^-) + /(^+)), (2.5) 
the ODEs can be written as 

Xk^u{xk), rhk = -{b- l)mk{u^{xk))- (2.6) 

In the Camassa-Holm case b — 2, this is a canonical Hamiltonian system gen- 
erated hy h — 5 X]j fc=i "^i "^fe Explicit formulas for the n-peakon 
solution of the Camassa-Holm equation were derived by Beals, Sattinger and 
Szmigielski [H [2] using inverse spectral methods, and the same thing for the 
Degasperis-Procesi equation was accomplished by Lundmark and Szmigielski 

[23 [ig. 

It requires some care to specify the exact sense in which the peakon solu- 
tions satisfy the PDE. The formulation (|2.2p suffers from the problem that the 
product mUx is ill-defined in the peakon case, since the quantity m = u — Uxx = 
2 X]r=i is ^ discrete measure, and it is multipHed by a function Ux which 

has jump discontinuities exactly at the points Xk where the Dirac deltas in the 
measure m are situated. To avoid this problem, one can instead rewrite (|2.ip as 

(1 - dl)ut + {b+l-dl) dx (i u") + dx ul) = 0. (2.7) 

Then a function u{x,t) is said to be a solution if 

• u{-,t) G W[^'j,^(R) for each fixed t, which means that u(-,i)^ and Ux{-,t)^ 
are locally integrable functions, and therefore define distributions of class 
D'CR) (i.e., continuous linear functional acting on compactly supported 
C°° test functions on the real fine R), 

• the time derivative ut{-,t), defined as the limit of a difference quotient, 
exists as a distribution in I?'(R) for all t, 

• equation (|2.7p . with dx taken to mean the usual distributional derivative, 
is satisfied for all t in the sense of distributions in V (R) . 

It is worth mentioning that functions in the space W^^^CR) are continuous, by 
the Sobolev embedding theorem. However, the term is absent from equation 
l|2.7p if 6 = 3, so in that particular case one requires only that u{-, t) e L^q^(R); 
this means that the Degasperis-Procesi can admit solutions u that are not 
continuous [61171126]. 

It is often appropriate to rewrite equation l|2.7p as a nonlocal evolution equa- 
tion for u by inverting the operator (1 — 9^), as was done in i9j for the 
Camassa-Holm equation. However, the distributional formulation used here is 
very convenient when working with peakon solution. 
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3 Novikov's equation 



The new integrable equation found by Vladimir Novikov is 

Ut - Uxxt + ^U^Ux = SuUxUxx + U^Uxxx, (3.1) 

which can be written as 

rrit + {mxU + 3mUx) u = 0, m = u — Uxx, (3-2) 
to highhght the similarity in form to the Degasperis-Procesi equation, or as 

(1 - dl)ut + (4 - dl) dx (i u^) + dx (i uul) + \ul^Q (3.3) 

in order to rigorously define weak solutions as above, except that here one re- 
quires that u(-,t) € Wj^f (R) for all i, so that and u^. are locally integrable 
and therefore define distributions in I?'(R); it then follows from Holder's in- 
equality with the conjugate indices 3 and 3/2 that um^ is locally integrable as 
well, and (|3.3p can thus be interpreted as a distributional equation. Since func- 
tions in Wj^'j?'(R) with p > 1 are automatically continuous, Novikov's equation 
is similar to the Camassa-Holm equation in that it only admits continuous dis- 
tributional solutions (as opposed to the Degasperis-Procesi equation, which has 
discontinuous solutions as well). 

Like the equations in the 6-family (|2.ip . Novikov's equation admits (in the 
weak sense just defined) multipeakon solutions of the form (|2.3p . but in this 
case the ODEs for the positions and momenta are 

Xk = u{xkf = e~l"^'=~^'l^ , 

■nik = -rUk u{xk) {ux{xk)) (3.4) 

^ rrij sgn(xfc - Xj) e-l^'=-^jl 

These equations were stated in [19j, where it was also shown that they constitute 
a Hamiltonian system Xk = {xk,h}, rhk — {mk,h}, generated by the same 
Hamiltonian h — ^ 12^j k=i i^j^k e~l^^^^'=l as the Camassa-Holm peakons, but 
with respect to a different, non-canonical, Poisson structure given by 

{xj,Xk} = sgn(xj - Xfc) (l - £"1^), 

{xj,mk} = rukE'^jk, (3.5) 
{mj,TOfc} = sgn{xj - Xk) mj-mkE^^., where Ejk = e"'^'"'''''. 

As will be shown below, i|3.4p is a Liouville integrable system (Theorem l4.7p : in 
fact, it is even explicitly solvable in terms of elementary functions (Theorem l9.ip . 
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4 Forward spectral problem 





In order to integrate the Novikov peakon ODEs, we are going to make use of 
the matrix Lax pair found by Hone and Wang [19], specified by the following 
matrix linear system: 

''O zm 1 

zm I ( V2 I , (4.1) 

^10 

UxZ^^ — u^mz 

— —Uj.z~^ — u'^mz I I ■02 I ■ (4.2) 

uz~^ 

(Compared with reference [19] we have added a constant multiple of the identity 
to the matrix on the right hand side of l|4.2p . and used z in place of A.) In 
the peakon case, when u = X]r=i ' the quantity m = u — u^x = 

2 X)r=i "^i is a discrete measure. We assume that xi < X2 < ■ ■ ■ < Xn 
(which at least remains true for a while if it is true for i = 0). These points 
divide the x axis into n + 1 intervals which we number from to n, so that 
the fcth interval runs from Xk to Xk+i, with the convention that xq = — oo and 
Xn+i — +00. Since m vanishes between the point masses, equation l|4.ip reduces 
to dxi^i — "03) dx'4'2 — and dxipa = 0i in each interval, so that in the fcth 
interval we have 

/Ak + z^ Ck e-^\ 
02 = 2zBk for Xk < X < Xk+i, (4.3) 

^03/ \Ake^-z^Cke~y 

where the factors containing z have been inserted for later convenience. These 
piecewise solutions are then glued together at the points Xk- The proper inter- 
pretation of (|4.ip at these points turns out to be that must be continuous, 
while ipi and -02 are allowed to have jump discontinuities; moreover, in the 
term zmip2, one should take tp2(x)6x^, to mean {ip2{xk)}Sxk- This point is fully 
explained in Appendix [b1 This leads to 



(4.4) 





1 + Xml 



where A = —z'^. 



We impose the boundary condition {Aq, Bq,Co) ~ (1,0,0), which is consistent 
with the time evolution given by l|4.2p for x < xi. Then all (Ak, Bk,Ck) are 
determined by successive application of the jump matrices Sk{X) as in l|4.4p . 
For X > Xn, equation (j4.2p implies that {A, B,C) := (A„, B„, C„) evolves as 

i = o, 5^^^, c^^-M±^^^AMA^ (4.5) 
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where M+ = X^fcLi "^fe 6^*° • Thus A is invariant. It is the (1,1) entry of the 
total jump matrix 

S{X) = Sn{X)...S2iX)Si{X), (4.6) 
and therefore it is a polynomial in A of degree n, 

Aix) = ± H.i^xr = (i - r) • ■ ■ - ^) ' (4-^) 

A:=0 ^ 1 / V 



where Hq = 1 (since 5(0) — /, the identity matrix), and where the other coeffi- 
cients Hi, ... , Hn are Poisson commuting constants of motion (see Theorems l4.2l 
and lO below). 

The first linear equation l|4.ip . together with the boundary conditions ex- 
pressed by the requirements that Bq = Cq = and An{X) = 0, is a spectral 
problem which has the zeros Ai, . . . , A„ of A{X) as its eigenvalues. (To be pre- 
cise, one should perhaps say that it is the corresponding values of z = ±\/—X 
that are the eigenvalues, but we will soon show that the Xk are positive, at least 
in the pure peakon case, and therefore more convenient to work with than the 
purely imaginary values of z; see (|4.19p below.) 

EHmination of V'l from l|4.ip gives dx'ip2 — zm^ps and {d^ — l)?/'3 = zrm(j2, 
and the boundary conditions above imply that (V'2,V'3) ~^ (0)0) as a: ^ — oo 
and i/^s — *■ as X ^ -l-oo. Using the Green's function — e~'^'/2 for the operator 
9^ — 1 with vanishing boundary conditions, we can rephrase the problem as a 
system of integral equations. 



^2ix) ^ z ip3{y)dm{y), 

r 1 

i;s.{x) = -z -e 1^ y\^2{y)dm{y), 



_ 2 

with integrals taken with respect to the discrete measure m = 2^"^^ miS^.. 
Here, there is again the problem of Dirac deltas multiplying a function ■!/'2 with 
jump discontinuities, and we take '4'2{x)5xt, to mean the average {tp2{xk))Sxkj 
in full agreement with the earlier definition of the singular term appearing in 
the spectral problem. Then 

{^2ixj)) = z 2^^p3{xk)mk +tp3{xj)mj J , 

\ k=l 



(4.9) 



k=l 

which can be written in block matrix notation as 



^If) - = i-lp T) m . ("«) 
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where 



P = diag(TOi, . . . ,m„), (4-11) 

E - iE,k)lk=i, where E^k = e-l"^-"'=l, 

r = , where T,fe = 1 + sgn(j - k). 

(In words, T is the lower triangular nxn matrix that has 1 on the main diagonal 
and 2 everywhere below it.) In terms of (^'2) alone, we have 

(*2) = -z2rp£;p(*2), (4.12) 

so the eigenvalues are given by = det(/ + z^TPEP) = det(/ - \TPEP), 
where of course / denotes the nxn identity matrix. Since the eigenvalues are 
the zeros of A{X), and since ^(0) = 1, it follows that 

A{\) = det(/ - \TPEP). (4.13) 

This gives us a fairly concrete representation of the constants of motion Hk, 
which by definition are the coefficients of A(X) (see l|4.7p ). and it can be made 
even more explicit thanks to the curious combinatorial result in Theorem 14.11 
We remind the reader that a fc x fc minor of an n x n matrix X is, by definition, 
the determinant of a submatrix Xij = {Xij)i^i_ j^^j whose rows and columns are 
selected among those of X by two index sets /, J C {1, . . . , n} with k elements 
each, and a principal minor is one for which I — J. Compare the result of 
the theorem with the well-known fact that the coefficient of s'' in det(J + sX) 
equals the sum of all principal k x k minors of X, regardless of whether X is 
symmetric or not. 

Theorem 4.1 ("The Canada Day Theorem"). Let the matrix T be defined as 
in (|4.1ip above. Then, for any symmetric nxn matrix X , the coefficient of 
in the polynomial det(/ + sTX) equals the sum of all k x k minors (principal 
and non-principal) of X . 

Proof. The proof is presented in Appendix El It relies on the Cauchy-Binet 
formula, Lindstrom's Lemma, and some rather intricate dependencies among 
the minors of X due to the symmetry of the matrix. □ 

Theorem 14.11 is named after the date when we started trying to prove it: 
July 1, 2008, Canada's national day. (It turned out that the proof was more 
difficult than we expected, so we didn't finish it until a few days later.) Summa- 
rizing the results so far, we now have the following description of the constants 
of motion: 

Theorem 4.2. The Novikov peakon ODEs p.4p admit n constants of motion 
Hi, ... , Hn, where Hk equals the sum of all k x k minors (principal and non- 
principal) of the nxn symmetric matrix PEP = {mjmkEjk)'^ ^^i. (See l|4.1ip 
for notation.) 
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Proof. This follows at once from (|47|) . I|4.13p . and Theorem |411 □ 

Example 4.3. The sum of all 1 x 1 minors of PEP is of course just the sum 
of all entries, 

n n 

Hi^ mjUikE-jk = Y m^mfee"!''^"'^'"!, (4.14) 

j,k—l j,k—l 

and the Hamiltonian of the peakon ODEs l|3.4p is h — \Hi. Higher order 
minors of PEP are easily computed using Lindstrom's Lemma, as explained in 
Section IA.3I in the appendix. In particular, the constant of motion of highest 
degree in the ruk is 

n— 1 n 

= det{PEP) = n (1 - ^h+i) n (4.15) 

Example 4.4. Written out in full, the constants of motion in the case n = 3 
are 

,2 



Hi ^ 777,2 + + 2777l7772-Bl2 + 2 777, 1777 3 £'13 + 27772 TTlg i?23 , 

i/2 - (1 - Ef^) ml 777^ + (1 - ^1^3) ml ml + {l~ Ej^) ml ml 

+ 2(^23 - -B12-B13) 777 ^ 7772 "^3 + 2(£;i2 - -B13-B23) ^1 7772 "^3 ' 
i/3 = (1 - ^?2)(1 - S23) "^2 "13- 



(4.16) 



From now on we mainly restrict ourselves to the pure peakon case when 
mk > for all k (no antipeakons) . Our first reason for this is that we can then 
use the positivity of Hi and Hn to show global existence of peakon solutions. 

Theorem 4.5. Let 

V ^ {xi < ■ ■ ■ < Xn, all mk > 0} (4.17) 

be the phase space for the Novikov peakon system p.4p in the pure peakon case. 
If the initial data are in V, then the solution (x(t),m(t)) exists for all f G R, 
and remains in V . 

Proof. Local existence in V is automatic in view of the smoothness of the ODEs 
there. By (|4.14p and l|4.15p . both Hi and Hn are strictly positive on V. Since 
777^ < Hi, all mk remain bounded from above. The positivity of Hn ensures 
that the 777^ are bounded away from zero, and that the positions remain ordered. 
The velocities Xk are all bounded by (X]™fc)^) hence < Xk < C for some 
constant C, and the positions Xk{t) are therefore finite for all t £ H. Since 
neither Xk nor 777fc can blow up in finite time, the solution exists globally in 
time. □ 

Remark 4.6. The peakon ODEs p.4p are invariant under the transformation 
(7771, . . . , 777„) 1-^ {—mi, . . . , —mn), so the analogous result holds also when all 
?T7fc are negative. 
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Theorem 4.7. The constants of motion Hi, ... , Hn of Theorem \4.2\ are func- 
tionally independent and commute with respect to the Poisson bracket (|3.5p . so 
the Novikov peakon system p.4p is Liouville integrahle on the phase space V . 

Proof. To prove functional independence, one should check that J :— dHi A 
dH2 A ... A dHn does not vanish on any open set in V. Since J is rational in 
the variables {nik, e^''}^^i, it vanishes identically if it vanishes on an open set, 
so it is sufHcient to show that J is not identically zero. To see this, note that 

Hk=ek{ml...,ml) + OiEpg), (4.18) 

where denotes the kth elementary symmetric function in n variables, and 
0{Epq) denotes terms involving exponentials of the positions Xj. It is well 
known that the first n elementary symmetric functions are independent (they 
provide a basis for symmetric functions of n variables [29]), and therefore the 
leading part of J (neglecting the 0{Epq) terms) does not vanish. Since the 
0{Epq) terms can be made arbitrarily small by taking the Xk far apart, we see 
that there is a region in V where J does not vanish, and we are done. 

To prove that the quantities Hk Poisson commute with respect to the bracket 
l|3.5p , it is convenient to adapt some arguments of Moser that he applied to the 
scattering of particles in the Toda lattice and the rational Calogero-Moser sys- 
tem [31] . The Poisson bracket of two constants of motion is itself a constant of 
motion, so {Hj,Hk} is independent of time. Consider now this bracket at a fixed 
point (x°, m'^) := (a;°, a;^, . . . , a;° , m°, TO2, . . . m^) G V which we consider as an 
initial condition for the peakon flow (x(t), m(i)), which exists globally in time 
by Theorem l4.5l Theorem l9.41 which will be proved later without using what we 
are proving here, shows that the peakons scatter as t ^ —00; more precisely, 
tends to 1/Xk, while the Xk move apart, so that the terms 0{Epq) tend to zero. 
(It should also be possible to prove these scattering properties directly from 
the peakon ODEs, along the lines of what was done for the Degasperis-Procesi 
equation in [28, Theorem 2.4], but we have not done that.) Thus, from (|4.18p . 
{iJ,,i/fe}(x",m") = {ff„i?fe}(x(t),m(t)) = \imt^^^{Hj,Hk}{-K{t),m{t)) = 
limt^_oo{ej, efe}(x(t), m(t)). Now the Poisson brackets of these symmetric 
functions are given by hnear combinations of the Poisson brackets {mj,mk} 
with coefficients dependent only on the amplitudes. However, from l|3.5p it 
is clear that {toj, TOfe}(x(i), m(t)) = 0{Epq) 0, from which it follows that 
{cj, efc}(x(t), m(t)) ^ as i — > —00, and hence {Hj, Hk}{x° = as re- 
quired. □ 

Remark 4.8. Since the vanishing of the Poisson bracket is a purely algebraic 
relation, the Hk Poisson commute at each point of i?^", not just in the region V. 

The Xk, which are defined as the zeros of ^(A), are the eigenvalues of the 
inverse of the matrix TPEP, since A{\) = det(/ — XT PEP). Another reason 
why we restrict our attention to the case with all mk > is that the matrix 
TPEP can then be shown to be oscillatory (see Section fA. 2 1 in the appendix), 
which implies that its eigenvalues are positive and simple. Consequently, the 
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Afc are also positive and simple, and for definiteness we will number them such 
that 

0<Ai<---<A„. (4.19) 

(For another proof that the spectrum is positive and simple, see Theorem 16. Ih 
Turning now to B = S{X)2i and C = S'(A)3i, we find from (gj]) and (|44l) 
that they are polynomials in A of degree n—1, with B{0) = M+ and C(0) — M^, 
where M+ = J2k=i e^* as before. This means that the two Weyl functions 

.(A).-lgl a... C(A)^-g| (4.20) 

are rational functions of order 0(1/A) as A — > cxd, having poles at the eigen- 
values Afc. Let bk and Ck denote the residues, 

fc=i fe=i " 

The time evolution of {A, B, C), given by (|4.5p . translates into 

^(A) = "^^^^"^"\ C(A) = -c.(0)^(A). (4.22) 

Comparing residues on both sides in (|4.22p gives 

; bk . , , bk ^ b,nbk , . 

bk = T-, Cfe = T- = 2^ T-^- (4.23) 

Afc /\fc — Am'^k 

This at once implies bk{t) = 6fe(0) e*/"^*" , and integrating Cfe(T) from r — — oo 
(assuming that Ck vanishes there) to t = i then gives 



c ~ bmbk 

_^ Am + Afc 



A purely algebraic proof of this relation between the Weyl functions, not relying 
on time dependence and integration, will be given below; see Theorem 16.11 We 
note the identities X^i Cfc/Afc = ^d]" ^fe/^fc)^ and X^i AfcCfc = ^d]" ^fe)^- 

The multipeakon solution is obtained as follows. The initial data Xfc(O), 
mfc(O) (for k = l,...,n) determine initial spectral data Afc(O), 6fe(0), which 
after time t have evolved to Afc(t) = Afc(O), bk{t) = 6fc(0)e*/'^'= (since the Afc are 
the zeros of the time-invariant polynomial ^(A), and since the bk satisfy (|4.23p ). 
Solving the inverse spectral problem for these spectral data at time t gives the 
solution Xk{t), mk{t). The remainder of the paper is devoted to this inverse 
spectral problem. 
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5 The dual cubic string 



Just like for the Camassa-Holm and Degasperis-Procesi equations, some terms 
in the Lax pair's spatial equation (equation l|4.ip in this case, repeated as (|5.ip 
below) can be removed by a change of both dependent and independent vari- 
ables. We refer to this as a Liouville transformation, since it is reminiscent of 
the transformation used for bringing a second-order Sturm-Liouville operator 
to a simple normal form. This simplification reveals an interesting connection 
between the Novikov equation and the Degasperis-Procesi equation, and allows 
us to solve the inverse spectral problem by making use of the tools developed 
in the study of the latter. 

Theorem 5.1. The spectral problem 



d 



zm{x) 









zm{x) 


-02 





/ ' 


Us 



^ I V2 j = zm{x) I ( ^2 I (5.1) 

on the real line a; G R, with boundary conditions 

^2{x) 0, as X ^ — oo, 
e^^Val^;) — > 0, as X — > —oo, (5.2) 
£"^^'3(2;) — > 0, as X +00, 
is equivalent (for z ^ 0), under the change of variables 

y = tanhx, 
4'i (y) = "01 (x) cosh a; — -03 (x) sinh x, 
02 (y) = zip2{x), 
03(y) = z'^ 03 (a;)/ cosh a; , 
g{y) = m{x) cosh'^a;, 

A = -z^, 

to the "dual cubic string" problem 




(5.3) 



(5.4) 

-A 

on the finite interval ~1 < y < I, with boundary conditions 

02(-l) =03(-l) = 03(1) - 0. (5.5) 

In the discrete case m — 2 to^ (^xj. , the relation between the measures m 

and g should be interpreted as 

n 

9iy) = ^9k3yk> yfc^tanhxfc, 5^ = 2 cosh Xfc . (5.6) 

fc=i 
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Proof. Straightforward computation using the chain rule and, for the discrete 

case, S^^ ^ ^{xk)Sy^. □ 

Remark 5.2. The cubic string equation, which plays a crucial role in the 
derivation of the Degasperis-Procesi multipeakon solution ^28j, is 

= -Ag^, (5.7) 

which can be written as a system by letting $ = (0i, 02, ^3) = (0, 4'yi4>yy)- 



d_ 

dy 




•J2 = 

h \-)^9{y) 0/ V03> 



(5.8) 



The duality between (|5.4p and (|5.8p manifests itself in the discrete case as an 
interchange of the roles of masses gu and distances 1^ = yk+i — yk] see Section [6l 
When the mass distribution is given by a continuous function g{y) > 0, the 
systems are instead related via the change of variables defined by 

j^-5(y) = ^, (5.9) 

dy g{y) 

where y and g{y) refer to the primal cubic string l|5.8p . and y and g{y) to the 
dual cubic string l|5.4p (or the other way around; the transformation l|5.9p is 
obviously symmetric in y and y, so that the dual of the dual is the original 
cubic string again). 

Remark 5.3. The concept of a dual string figures prominently in the work of 
Krein on the ordinary string equation dy(f> = —Xgcj) (as opposed to the cubic 
string). For a comprehensive account of Krein's theory, see |12j . 



Remark 5.4. As a motivation for the transformation l|5.3p . we note that one can 
eliminate from (|5.ip . which gives = zmili^, {d^ — 1)1^3 — zm^2- From 
the study of Camassa-Holm peakons |i2j it is known that the transformation 
y = tanha;, 4'{y) = ipix)/ cosh a; takes the expression (3^ — 1) t/i to a multiple 
of (j)yy , so it is not far-fetched to try something similar on ip3 while leaving ip2 
essentially unchanged. 

From now on we concentrate on the discrete case. The Liouville transfor- 
mation maps the piecewise defined ("01! '/'2, V'a) given by (|4.3p to 



"1 



Ak{X)-XCkiX) 



h\ = \ -2XBk{X) for yk<y<yk+i. (5.10) 



-XAk{X){l+y)-X^Ck{X)il~y) 



The initial values {Aq, Bq, Cq) = (1,0, 0) thus correspond to $(-1; A) = (1, 0, 0)*, 
where $(y; A) = (0i, 02, ^a) , and at the right endpoint y = 1 we have 

'An{X)-XCr.{X)' 

$(1;A)=( -2AB„(A) |. (5.11) 
-2AA„(A) 
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In particular, the condition An{X) = defining the spectrum corresponds to 
(^3(1; A) = 0, except that the latter condition gives an additional eigenvalue 
Ao = which is only present on the finite interval. (This is not a contradiction, 
since the Liouville transformation from the line to the interval is not invertible 
when z = -A^ = 0.) 

The component 03 is continuous and piecewise linear, while (fii and 02 are 
piecewise constant with jumps at the points yk where the measure g is supported. 
More precisely, at point mass number k we have 



(5.12) 



'1 


9k 






1 


9k 











1 


1 


I 





1 1 




-A/fe 






4)i{yt) - Mvk ) gkiMvk)), 

02(2/fe ) - hiVk ) = 9k hiVk), 
and in interval number fc, with length Ik ~ yk+i — Vk, 

Mvk+i) - Mvt) = hdyUvi) = --xikMyi)- (S-is) 

In terms of the vector <& these relations take the form 

A 9k yi\ 

'^(yt) - I 1 9k \ $(yfe ), (5.14) 

and 

*(2/fe+i)= (^0^^ 0]^y+), (5.15) 

respectively. If we introduce the notation 
/I 0\ 

G{x,X)=\ 1 0|, L(a;) = I 1 x \ , (5.16) 

\-Xx 

it follows immediately that 

<i>(l;A) = G(;„,A) L(.9„) ...G{h,X) Lig^) G(/i,A) L{g,) G(/o,A)(o) . 

(5.17) 

We define the Weyl functions W and Z of the dual cubic string to be 

urtw '?^2(i;A) A) 

^ - , Z{X) = - 5.18) 

03(1; A) 03(1; A) 

It is clear from (|5.1ip that they are related to the Weyl functions lu and C 
previously defined on the real line (see l|4.20p ) as follows: 





X 


1.^ 





1 


X 








1 



k 



(5.19) 

_ A„(A)-AG„(A) _ 1 1 , Cfc 

^^^^ ~ 2AA„(A) ~ 2A + ^^^^ - 2A + ^ A^' 



13 



6 Relation to the Neumann-like cubic string 



Kohlenberg, Lundmark and Szmigielski |23j studied the discrete cubic string 
with Neumann-like boundary conditions. We will briefly recall some results 
from that paper, with notation and sign conventions slightly altered to suit our 
needs here. The spectral problem in question is 

'l^yyyiy) = ~^9{y)4>{y) for y e R, 

(j)y{-00) = (j)yy{~00) = 0, 4>yy{oo) = 0, 

where g = X]fc=o 9k ^yk ^ discrete measure with n + 1 point masses g^, . . . ,gn 
at positions j/o < yi < • ■ ■ < Vn', between these points are n finite intervals of 
length h, . . . Jn (where h = Vk — Vk-i)- Since (j)yyy = away from the point 
masses, the boundary conditions can equally well be written as 

MVo ) = 'l^yyiVo ) "^^^ 'l^yy iVn) ^ 0. 
Using the normahzation 0(— oo) = 1 (or (t>{y^) = 1) and the notation $ = 

(0, 4)y,4>yyY, 0116 fittdS 



<&( 



>(2/+;A) -G(5„,A) L(/„) . . . G(g2, A) Lih) G(5i,A) L{h) G{go,\) (o) , 

(6.2) 

with matrices G and L as in (|5.16p above. Under the assumption that all gk > 0, 
the zeros of (j^yyiyt.] A), which constitute the spectrum, are 

= Ao < Ai < • • • < An, 

and the Weyl functions are 



(l)vy{yn;X) 7A ^A-Afc ^ 



(6.3) 



with all bk > 0. They satisfy the identity 

Z{X) + Z{-X) + W{X)Wi-X)=Q, (6.4) 
from which it follows, by taking the residue at A = Afc, that 

Ebmbk _x 
ir"TAr- (^-^^ 

m— 1 

Thus Z(X) is uniquely determined by the function W{X) and the constant 7. 

Now note that (|6.2p is exactly the same kind of relation as (|5.17p , except that 
the roles of ^fe and 1^ are interchanged, and the right endpoint is called y ^ 
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instead of y = 1. The definitions of the Weyl functions (|6.3p also correspond 
perfectly to the Weyl functions l|5.18p for the dual cubic string. Therefore, all 
the results above are also true in the setting of the dual cubic string. The 
assumption that the n distances Ik and the n + 1 point masses gk are all positive 
for the Neumann cubic string corresponds of course to the requirement that the 
n point masses gk and the n+1 distances Ik are positive for the dual cubic 
string. The constant 7 = J2k=o 9k in the term I/7A in (|6.3p corresponds to the 
constant 2 in the term 1/2A in l|5.19p . since X)fe=o ^'^ = 2 is the length of the 
interval — 1 < y < 1. In summary: 

Theorem 6.1. Assume that all point masses gk are positive. Then the discrete 
dual cubic string of Theorem \ 5.1\ has nonnegative and simple spectrum, with 
eigenvalues = Aq < Ai < • • • < A„, and its Weyl functions l|5.18p have 
positive residues and satisfy (|6.4p and l|6.5p . In particular, the second Weyl 
function Z{\) is uniquely determined by the first Weyl function W{X). 

7 Inverse spectral problem 

The inverse spectral problem for the discrete dual cubic string consists in re- 
covering the positions and masses {Vk, gk}k=i given the spectral data consisting 
of eigenvalues and residues {Afc,6fc}^^^ (or, equivalently, given the first Weyl 
function W{Xj). The corresponding problem for the Neumann-like cubic string 
was solved in [23|, and we need only translate the results, as in Section [6l See 
also |28] for more information about inverse problems of this kind and [3J for 
the underlying theory of Cauchy biorthogonal polynomials. 

To begin with, we state the result in terms of the bimoment determinants 
T>m''^ and I?^ defined below. Things will become more explicit in the next 
section fCorollarv l8.4p . where the determinants are expressed directly in terms 
of the Afc and bk- 

Definition 7.1. Suppose ^ is a measure on R+ (the positive part of the real 
line) such that its moments. 



and its bimoments with respect to the Cauchy kernel K{x, y) ~ \/{x + y), 



are finite. For to > 1, let T>m denote the determinant of the mx m bimoment 
matrix which starts with lab in the upper left corner: 




(7.1) 




(7.2) 



,{ab) 





,(ba) 



(7.3) 



m 




a-\-7n— l,b-\-7n— 1 
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Let Vi"'''^ = 1, and V, 



(ab) 



for TO < 0. 



Similarly, for m > 2, let denote the m x m determinant 



f3o 


ho 


hi ■ 


h,m-2 


Pi 


ho 


hi . 


h,m-2 


P2 


ho 


/31 ■ 


h,m~2 



(7.4) 



Pm — l lynO hnl • ■ ■ hn^m — 2 

and define !)[ = i3q and = for m < 1. 

Theorem 7.2. Given constants < Ai < • • • < A„ and 5i, . . . , 5„ > 0, define 
the spectral measure 



and let hb be its bimoments, 



hb — 



i=l J=l 



" ' b,b 



Ai + Xj 



(7.5) 



(7.6) 



Then the unique discrete dual cubic string (with positive masses gu) having the 
Weyl function 



f-fA-Afe J X-K 



is given by 



Vk' 



k=l 



(00) _ 1^(11) 

7^(00) , i-n(ii) 



V 



gk' = 2 



T-,(00) , 1^(11) 
^k "I" 2^k-l 



(7.7) 



where k' = n+1 — k for fc = 0, . . . , n + 1. The distances between the masses are 
given by 



h'-i = Vk' - Uk' 



(11) 
fe 



(7.8) 



Proof. For < fc < n, let a(^'^+^^(A) be the product of the first 2k + 1 factors 
in (ICTtI) . 



a(2^-+i)(A) =G(?„,A) L(g„) G(;„_i,A) L(.9„_i) . , 



G{lk',X) L{gk') G{lk'-i,X), (7.9) 
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where k' = n + 1 — k. By Lemma 4.1 and Theorem 4.2 in [23| . the entries in the 
first column of a — a(^'^+^^(A), 



'an' 

021 



satisfy what in [23] was called a "Type I" approximation problem. This means 
that (-P(A), -P(A), (5(A)) are polynomials in A of degree k, k, k + 1, respectively, 
satisfying the normalization conditions 

P(0) = 1, P(0) = 0, Q(0) = 0, 

the approximation conditions 

Q(A)M^(A) +P(A) = 0(1), Q(A)Z(A) +P(A) = 0(A-^), as A c5o, 

and the symmetry condition 



Q(A)Z(-A) - P(A)W^(-A) - P(A) = O(A-^-i), 



as A — > oo. 



According to Theorem 4.15 in [23], this determines {P,P,Q) uniquely; in par- 
ticular, the coefficients of a^^i^^\X) — Q{\) — J2i=i Qi^^ given by the 
nonsingular linear system 



ho 
ho 



From l|7.9p one finds that 

(2fc+l) 



loi ■ 


■ hk\ 


( qi\ 




hi ■ 


■ hk 







hi ■ 


hk 


93 = 





hi ■ 


IkkJ 


\9fc+l/ 





(7.10) 



•^31 



(A) = {-\){ln + ln-1 + ■■■+ h'-l) + 



-A) 



2 2 
k+l ( 9n 9n~l 

2 2 



9l 



-11 1 



(7.11) 



and the lowest and highest coefficients are then extracted from l|7.10p using 
Cramer's rule: 



V 



(11) 



-91 



7-,(00) l7-,(ll) 

'^k+1 + 2-^fc 



^3 l-2/fc'-i, 

]=k'-l 



(,-lf+\k+l = 



V 



(10) 



r)(00) 1^(11) 
■^fc+1 + 2-^fe 



(7.12) 



h'-i- 



The first equation gives a formula for yk'-i right away, and of course also for yk' 
(with 1 < k < n + 1) after renumbering. This formula l|7.7p for yk' holds also 
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for k = 0, since it gives yo' = Un+i = +1 because of the way I'm''-* is defined 
for m < 0. (That it indeed gives y(n+i)' — Uo = when fc = n + 1 is not 
as obvious; this depends on I?^*^! being zero when the measure ^ is supported 
on only n points. See [23l Appendix B].) Subtraction gives a formula for Ik'-i 
which simplifies to (|7.8p with the help of "Lewis Carroll's identity" [Ml Prop. 10] 
appHed to the determinant 2?^°^] : 

= - ^ . (7.13) 

Finally, the second formula in (|7.12p , divided by the corresponding formula with 
k replaced by fc — 1, gives an expression for ^ gl, Ik'-i from which one obtains 




^(lO)p(lO) ■ 

The formula for gk' presented in l|7.7p now follows from the identity (2?^)^ = 
2pUo)p(^io) ^j^g positivity of I?^, which are immediate consequences of (|8.6p 
below. (The determinant identity can also be proved directly by expanding V^. 
along the first column, squaring, and using Pif3j — li+ij + lij+i.) □ 

Remark 7.3. We take this opportunity to correct a couple of mistakes in [23] : 
the formula in Corollary 4.17 should read [Q3k+2] = {—^)''^^T>k/Ak+i, and 

consequently it should be Wn-fc = 77-; — ;— in (4.54). 

2Ak+iAk 

8 Evaluation of bimoment determinants 

The aim of this section is just to state some formulas for the bimoment determi- 
nants Vm''^ and V^, taken from |28[ Lemma 4.10] and |23[ Appendix B]. Quite 
a lot of notation is needed. 

Definition 8.1. For fc > 1, let 

1 f A(x)2 dfi'^ix) 



tk 



kl J^k T{x) xiX2---Xk 



1 f A a; 2 
kl J^k T{x) 

fc! J^k T{x) 



where 



A(a;) = A(xi, ...,Xk)= Y[i^i - ^j)' 
V{x) = r(a;i, . . . ,Xfe) = Y[.^^i + ^j)- 



i<3 
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(When fc = or 1, let A(a;) = r(a;) = 1.) Also let to = uq = vq = 1, and 

tk = Uk = Vk = for fc < 0. 

When fj, = J2k=i Sx^, the integrals tk, Uk, Vk reduce to the sums Tk, Uk, 
Vk below. 

Definition 8.2. For fc > 0, let ('^^"^) denote the set of fc-element subsets / = 
{ii < ■ ■ ■ < ik} of the integer interval [1, n] = {1, . . . , n}. For / G ('^fc"') ; let 

A/ = A(A,,,...,A,J, T/ = r(A,,,...,A,J, (8.3) 

with the special cases A0 = = A^jj = F^jj. = 1. Furthermore, let 

iei iei 
A? 

with A0 = 60 = 1. Using the abbreviation = —, let 

1 1 



A/ 



and 



Wk 
Zk 



Uk = 


E 






Uk 


Vk-1 


Uk+1 


Vk 


Tk 


Uk-1 


Tk+i 


Uk 



= UkVk — Uk+iVk-1, 
= TkUk - Tk+iUk-i- 



(8.5) 



(To be expHcit, Ua = Vo = To = 1, and [/fe = Vfc = Tfc = for fc < or fc > n.) 

We can now finally state the promised formulas for the bimoment determi- 
nants. 

Lemma 8.3. For all m, 



tjji Ufn — l 




^m— 1 




.. 2)(ii) = . 





2)(io) 



(Um) 



v = 



2m 



(8.6) 



2m ' m 2m— 1 

/n ifte discrete case when ji = bk S\^, , this reduces to 



fe=i 



~(00) _ ^rn -r,(ll) ^ '^m r)(10) = i^'mY 
m 2m ' m 2m ' m 2m 



V' 



Ura Um — 1 



•7) 
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Corollary 8.4. The solution to the inverse spectral problem for the discrete 
dual cubic string (Theorem \7.2\) can be expressed as 



Zk - Wk-i Zk + Wk-i ,^ 

h'-, = yu' - yu'-, - ^z, + W,^,){Z,^, + W,y 

The expression Wk can be evaluated explicitly in terms of and bk, al- 
though the formula is somewhat involved [28t Lemma 2.20]: 



m=lT^([-iM\ y \ CUD=J 

^U-mJ |C| = |i3|=m 

J^C'l'Z^) min(C)<min(_D) 

/nj=0 



(8.10) 



where Aj j = (Aj — Aj)^. The corresponding formula for Zk is obtained 

iei.jeJ 

by replacing 6^ with bi/Xi everywhere. 



9 The multipeakon solution 

In order to obtain the solution to the inverse spectral problem on the real line, 
which provides the multipeakon solution, we merely have to map the formulas 
for the interval (Corollarv l8.4p back to the Hne via the Liouville transformation 

We remind the reader that in this paper we primarily study the pure peakon 
case where it is assumed that all ruk > and also that xi < • ■ ■ < x„. This 
assumption guarantees that the solutions are globally defined in time (Theo- 
rem l4.5p and, regarding the spectral data, that all bk > and < Ai < • • • < A„ 
(Theorem lG.ip . Details regarding mixed peakon- antipeakon solutions are left for 
future research, but we point out that since the velocity Xk — u{xk)'^ is always 
nonnegative, Novikov antipeakons move to the right just like peakons (unlike 
the 6-family l|2.ip , where pure peakons move to the right and antipeakons to the 
left, if they are sufficiently far apart). Nevertheless, peakons and antipeakons 
may collide after finite time also for the Novikov equation, causing division by 
zero in the solution formula for in l|9.ip below, and this breakdown leads 
to the usual subtle questions regarding continuation of the solution beyond the 
collision. 
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Theorem 9.1. In the notation of Section[^ the n-peakon solution of Novikov's 
equation is given by 



1 , Zk 

— m . 

2 Wk-i 



(9.1) 



UkUk-i 

where k' = n + I — k for k = I, . . . ,n, and where the time evolution is given by 

bk{t)^bk{0)e'/^''. (9.2) 
Proof. The inverse of the coordinate transformation (|5.6p is 

1 



Xk 



-mi- 

2 1 



Vk 
Vk 



ruk 



9k^^-yl 



which upon inserting l|8.8p gives (|9.ip at once. The evolution of bk comes from 
equation (|423l) . □ 

Example 9.2. The two-peakon solution is 

(Ai - x^Y 



1, Z-2 

xi = — m 

^ 2 Wi 



1 



62 /j2 

(Ai+A2)2AiA2 ' ' 
X,bl + X,bl + 4^b,b, 



■In- 



X2 



2 Wo 



■In 



1 
Ai 



Ao 



Ai + A2 
4 

■bib, 



Ai + A: 



mi = 



U2U1 



(Ai-A2)H26i 



(Ai + A2)2AiA2 



Ai 6? + A2 52 + 1^6162 
Ai + A2 



1/2 



(Ai- A2)^&i62 



Xibl 



X^bl 



Ai + A2 

4A1A2 
Ai + A2 



(foi + 62) 

1/2 



(9.3) 



&162 



VA1A2 (fei + 62) 



m2 



6? 



V^iW^o \Xi A2 Ai + A: 



■ 6165 



1/2 



where the simpler of the two expressions for mi is obtained under the assumption 
that all spectral data are positive, and therefore only can be trusted in the pure 
peakon case. This way of writing the solution is simpler and more explicit 
than that found in [19]. In order to translate i|9.3p to the notation used there, 
write {qk,Pk) instead of {xk,rak), Ck instead of 1/Afc, and to instead of {X^^ — 
X^^)-^ In then tanhT = (6i-62)/(&i+^>2) and cosh"^ T = 46l62/(6l+62)^ 
where T — ^(ci — C2){t — to). 
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Example 9.3. The three-peakon solution is 

^^=2^"W^' "^=2^"^!' "^=2^"in^ 



A1A2 A1A3 A2A3 

4«'i3«'23 .,,2_L 4^-12*23 ,,2i,_l 4 *i2«'i3 ,2,, 

7 _ ^123 1,2. 2. 2 
A1A2A3 



and 



^ (Ai-A2)^ (Ai-A3f (A2-A3)2 
*12 = — ; — — ; , ^13 = — r — ; , ^^23 



A + A2 Ai + A3 A2 + A3 

(Ai-A2)2(Ai-A3)2(A2-A3)2 



123 



(Ai + A2)(Ai +A3)(A2 + A3) 



(9.4) 



VZ3W2 VZ2W1 VZiWo 

where Uo = Wo = 1, 

Ui=bi+b2 + h, 

U2 = *12 hb2 + *13 bibs + *23 &2&3, (9.5) 
Us = *123&1&2&3, 

= Ai 6? + A2 bl + A3 6^ 

, 4A1A2 , , , 4A1A3 , , , 4A2A3 , , 

■O1O2 + -; — ^0l03 + T —^0203, 



Ai + A2 Ai + A3 A2 + A3 

W2 = *?2 A1A2 bftl + A1A3 bjbl + ^^3 A2A3 blbj ^9 6) 

4^-13*23 A1A2A3 , , , 2 , 4*i2*23AiA2A3 3, 

+ A,+A2 '^'^'^ + A,+A3 '''''' 

4^'i2«'l3 A1A2A3 ,2u u 
+ A2+A3 '^'^'^^ 

Z _ ^1 I I I 4 ^ ^ ^ 4 ^ ^ ^ 4 ^ ^ 
Ai A2 A3 Ai + A2 Ai + A3 A2 + A3 



(9.7) 



(9.8) 



Theorem 9.4 (Asymptotics). Let the eigenvalues he numbered so thatO < Ai < 
• • • < A„ . Then 

t 1 " fA — A 1^ 

Xk{t) ^ — + logbkiO)- -In \k+ Yl 1^7T^~r~rTr' as -00, 
Afe 2 (Ai + Afe)Ai 

(9.9) 



Xk'{t) ~ — + log6fc(0)- - In Afc + V In ' \ , as i ^ +00, 
Afc ^ f— f (Aj + AfcJAj 
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where k' = n + I — k. Moreover, 



lim nikit) = = lim mk'{t). (9.10) 



In words: asymptotically as t ~* ±00, the kth fastest peakon has velocity 1/Xk 
and amplitude l/^/Xk- 



Proof. This is just a matter of identifying the dominant terms; bi(t) — 61 (0) e*/'^^ 
grows much faster ast~+ +00 than 62 (t) , which in turn grows much faster than 
b3{t), etc., and as t ^ —00 it is the other way around. Thus, for example, 
Uk ~ ^ 12. . .kbib2 ■■ - bk as t +00. A similar analysis of Wk and Zk leads 
quickly to the stated formulas. □ 

The only difference compared to the Xk asymptotics for Degasperis-Procesi 
peakons [281 Theorem 2.25] is that l|9.9p contains an additional term — ^InAfc. 
Since this term cancels in the subtraction, the phase shifts for Novikov peakons 
are exactly the same as for Degasperis-Procesi peakons j28] Theorem 2.26]: 

lim (xk'{t) - T~ ) ~ 1™ (xkit) - 

t^oo \^ \k J \ Afe 

E^i (Ai - Afe)^ >^ (Ai - Afc)^ /r,11^ 
._/°^(a7Ta;)a--.5/°''(a7Ta^ ^'-''^ 



A Combinatorial results 

This appendix contains some material related to the combinatorial structure of 
the constants of motion Hi, ... , Hn of the Novikov peakon ODEs; see SectionlH 
and in particular Theorem 14.21 Recall that 

^(A) = 1 - Affi + • • • + (-A)"ff„ det(/ - XTPEP), 

where / is the n x n identity matrix, and T, E, P are n x n matrices defined 
by Tjk = 1 + sgn(j - k), Ejk = e-I^J"^"!, and P = diag(mi, . . . ,m„). The 
first thing to prove is that the matrix TPEP is oscillatory if all ruk > 0, which 
shows that the zeros of ^(A) are positive and simple. Then we show how to easily 
compute the minors of PEP, and finally we prove the "Canada Day Theorem" 
(Theorem 14. ip which implies that Hk equals the sum of all kxk minors of PEP. 



A.l Preliminaries 

In this section we have collected some facts about total positivity \21\ fTSj fT3] 
that will be used below. 

Definition A.l. If X is a matrix and / and J are index sets, the submatrix 
{Xij)i^ij^j will be denoted by Xjj (or sometimes Xjj). The set of fc-element 
subsets of the integer interval [1, n] = {1, 2, . . . , n} will be denoted ('\"'), and 
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elements of such a subset / will always be assumed to be numbered in ascending 
order ii < ■ • ■ < ik- 

Definition A. 2. A square matrix is said to be totally positive if all its minors 
of all orders are positive. It is called totally nonnegative if all its minors are 
nonnegative. A matrix is oscillatory if it is totally nonnegative and some power 
of it is totally positive. 

Theorem A. 3. All eigenvalues of a totally positive matrix are positive and of 
algebraic multiplicity one, and likewise for oscillatory matrices. All eigenval- 
ues of a totally nonnegative matrix are nonnegative, hut in general of arbitrary 
multiplicity. 

Theorem A. 4. The product of an oscillatory matrix and a nonsingular totally 
nonnegative matrix is oscillatory. 

Definition A. 5. A planar network (r,w) of order n is an acyclic planar directed 
graph r with arrows going from left to right, with n sources (vertices with 
outgoing arrows only) on the left side, and with n sinks (vertices with incoming 
arrows only) on the right side. The sources and sinks are numbered 1 to n, 
from bottom to top, say. All other vertices have at least one arrow coming in 
and at least one arrow going out. Each edge e of the graph F is assigned a 
scalar weight i^{e). The weight of a directed path in F is the product of all the 
weights of the edges of that path. The weighted path matrix f2(F, lo) is the n x n 
matrix whose entry fiy is the sum of the weights of the possible paths 

from source i to sink j. 

The following theorem was discovered by Lindstrom [25] and made famous 
by Gessel and Viennot [16j. A similar theorem also appeared earlier in the work 
of Karlin and McGregor on birth and death processes p2]. 

Theorem A. 6 (Lindstrom's Lemma). Let I and J be subsets o/{l, . . . , n} with 
the same cardinality. The minor detfi/,/ of the weighted path matrix VL{T,u)) 
of a planar network is equal to the sum of the weights of all possible families of 
nonintersecting paths (i.e., paths having no vertices in common) connecting the 
sources labelled by I to the sinks labelled by J. (The weight of a family of paths 
is defined as the product of the weights of the individual paths.) 

Corollary A. 7. // all weights ^(e) are nonnegative, then the weighted path 
matrix is totally nonnegative. 

Remark A. 8. Beware that having positive weights does not in general imply 
total positivity of the path matrix fi, since some minors det flu may be zero 
due to absence of nonintersecting path families from / to J, in which case O is 
only totally nonnegative. 

A.2 Proof that TPEP is oscillatory 

The matrix T is the path matrix of the planar network whose structure is 
illustrated below for the case n = A (with all edges, and therefore all paths and 
families of paths, having unit weight): 
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© — — 
® — — 



Indeed, there is clearly one path from source i to sink j ii i — j, two paths if 
i > j, and none if z < j, and this agrees with 



1, 


i =j, 


2, 


i > j, 


0, 


i < j. 



Similarly one can check that the matrix PEP is the weighted path matrix of 
the planar network illustrated below for the case n = 5 (we are assuming that 
xi < ■ ■ ■ < Xn, so that E12E23 = 6^=1-^25^2-2:3 ^ g^p ). 

©1715 1 - m5 ^ 
^o- y .@ 

@^^^^ 5l>o i-^l ^@ 

©7712 -£-23^ 1 — EI2 /^^23 7712 
^o- .© 

©^I^^ ^© 

By Corollary IA.7[ both T and PEP are totally nonnegative (if all Wfc > 0). 
Furthermore, {PEP)^ is the weighted path matrix of the planar network ob- 
tain by connecting N copies of the network for PEP in series, and if N is large 
enough, there is clearly enough wiggle room in this network to find a nonin- 
tersecting path family from any source set / to any sink set J with |/| = |J|. 
Thus {PEP)^ is totally positive for sufficiently large N; in other words, PEP 
is oscillatory. (Another way to see this is to use a criterion |15[ Chapter II, 
Theorem 10] which says that a totally nonnegative matrix X is oscillatory if 
and only if it is nonsingular and Xij > for |« — j| = 1.) Since T is nonsingular. 
Theorem IA.4I impfies that TPEP is oscifiatory, which was the first thing we 
wanted to prove. 
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A.3 Minors of PEP 



Having a planar network for PEP makes it easy to compute its minors using 
Lindstrom's Lemma. 

Example A. 9. Consider the constant of motion in the case n = 6. 

For sources / = {1, 2, 3} and sinks J — {1, 2, 3} there is only one family of 
nonintersecting paths, namely the paths going straight across. The weights of 
these paths are mimi, TO2(1 — £'i2)"^2 and 7713(1 — £"13)7713, and the total weight 
of that family is therefore (1 — £i2)(l — £'23) rn\rn^m\^ which will be the first 
term in H^. 

A similar term results whenever I — J. For instance, when / = J = {1, 2, 4} 
the paths starting at sources 1 and 2 must go straight across, while the path 
from source 4 to to sink 4 can go straight across, or down to line 3 and up again. 
The contributions from these two possible nonintersecting path families add up 
to 

'2 \„, /I 772 \™ I ™ IP /'I 1^2 



777 i77li • 7712(1 — £12)^2 ■ (^7714(1 — £34)7714 + 7774£34(1 — £'23)-^34™4 

= {1-Ef2)il-El^)mlmlml 

From / = {1, 2, 3} to J = {1, 2, 4} there is one nonintersecting path family, 
and there is another one with the same weight from / = {1, 2, 4} to J = {1, 2, 3}; 
the two add up to the term 2(1 — £i2)(l — £|3)£24 777^771377737714. 

Continuing hke this, one finds that the types of terms that appear in H3 are 

iJg = (1 - £?2)(1 - -BI3) mlmlml + ... 

+ 2(1 - E^2)i'^ - £^23)^34 mlmlm^mi + ... 

+ 4(1 - £i2)(l - £|4)^23£^45 m\m2mzmim^ + . . . 

+ 8 (1 - £23)(1 ^ Elr,)Ei2EuE5(j 171117121713171411151125. 

The last term comes from the 8 possible nonintersecting path families from 

I = {ii,i2,i3} to J = {ji,j2,j3} where (ii, ji) = (1,2) or (2,1), (12,^2) = (3,4) 
or (4,3), and (13, J3) = (5,6) or (6,5). 

Remark A. 10. Alternatively, the 777^ can be factored out from any minor of 
PEP, leaving the corresponding minor of E, which can be computed using a 
result from Gantmacher and Krein [TSt Section II. 3. 5], since the matrix E is 
what they call a single-pair matrix. This means a symmetric n x n matrix X 
with entries 

X.,^h^- (A.2) 



The kxk minors of a single-pair matrix are given by the following rule: det Xjj 

k 



0, unless ('^l"') satisfy the condition 



(iiji) < (i2, J2) < • • • < {ik,jk), (A.3) 
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where the notation means that both numbers in one pair must be less than both 
numbers in the following pair; in this case, 



det Xij = Vai 



X/3i Xa2 




Xp2 


Xaa 






#1 V'as 




^/32 









X/3fc 



where 



,l3ra) = (min(i™, j„),max(i„, 



(A.4) 
(A.5) 



In the case of E we have -0* = e"^' and Xi = e (assuming as usual that 
xi < ■ ■ ■ < Xn), and l|A.4p becomes 

det Ei J = {l-El^Jil-El^J . . . {l-El_^^jE^,f3,E^,0, . ..E^,p,. (A.6) 
A.4 Proof of the "Canada Day Theorem" 

The result to be proved (Theorem l4.ip is that for any symmetric nxn matrix X, 
the coefficient of s'^ in the polynomial det(/ + s TX) equals the sum of all fc x /c 
minors of X: 



det(/ + sTX) = 1 + ^ I ^ dctX/j I s- 



(A.7) 



k=l 



We start from the elementary fact that for any matrix Y , the coefficients in its 
characteristic polynomial are given by the sums of the principal minors, 



det(/ + s r) = 1 + ^ Yjj s\ 

Applying this to y = TX and computing the minors of TX using the Cauchy- 
Binet formula [H Ch. I, § 2] 



det{TX)AB = ^-4/ dot XjB, for A, Be {^\"^) 



{A.S 



we find that 



det(/ + s TX) = 1 + XI I Tj/ det Xjj | s* 



k=l 



Comparing this to (jA.Tp . it is clear that what we need to show is that, for any k, 
Y Y det Tj/ det X/,7 = Y Y detX/j. (A.9) 



The first thing to do is calculate the minors det Tj/. 
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Definition A. 11. Given I,J £ ('^i"'); the set / is said to interlace with the 
set J, denoted I < J, if 

h<ji<i2<j2<---<ik<jk- (A. 10) 

If all the inequalities are strict, then / is said to strictly interlace with J, in which 
case we write I < J. If / < J, then /' and J' denote the strictly interlacing 
subsets (possibly empty) 

/' = /\(/nJ), J' = J\(/nJ), (A. 11) 

whose cardinality (possibly zero) will be denoted by 

p{I,J) = \I'\^\J'\. (A.12) 
Lemma A.12. For I,JE ('\"'); the corresponding k x k minor of T is 

detT,,^ ^ ^f'^'' (A.13) 

0, otherwise- 



Proof. We will use Lindstrom's Lemma (Theorem IA.6P on the planar network 
for T given in Section IA.2I above; the minor det Tjj equals the total number 
of families of nonintersecting paths connecting the source nodes (on the left) 
indexed by J to the sink nodes (on the right) indexed by /. 

The proof proceeds by induction on the size n of T. The claim is trivially 
true for rt = 1. Consider an arbitrary n > 1, and suppose the claim is true for 
size n — 1. If neither / nor J contain n, the claim follows immediately from the 
induction hypothesis, and likewise if / and J both contain n, because there is 
only one path connecting source n to sink n. If / contains n but J does not, 
then det Tji = because there are no paths going upward; this agrees with the 
claim, since in this case / does not interlace with J. 

The only remaining case is therefore J = JiU{n}, / = IiU{ik}, with ik < n. 
But then 

2, if jfc-i < ik, 
det Tj/ = detTji/i x S 1, if jfc-i = 

0, if jfe-i > ik, 

depending on whether the path connecting source n with sink ik has to cross 
the jk-i level; if it does not, there are two available paths, if it does, there is 
only one available path provided jk-i = in, otherwise the path intersects the 
path coming from source jk-i- In the last instance, / does not interlace with J, 
while in the other two / < J if and only if Ii < Ji, thus proving the claim. □ 

According to this lemma, the structure of l|A.9p (which is what we want to 
prove) is 

J2 2P^'^''> detXjj= detX^s, (A.14) 

/,je(i^r') ^.se(i^-i) 

KJ 
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and we must show that those terms det Xjj that occur more than once on the 
left-hand side exactly compensate for those that are absent. This will follow 
from another technical lemma: 

Lemma A. 13 (Relations between kxk minors of a symmetric matrix). Suppose 
I,J £ and I<J. Then, for any symmetric n x n matrix X , 

det Xab = 2"'^^'-^^ det Xjj. (A. 15) 

AnB=inJ 

Before proving Lemma IA.131 we will use it to finish the proof of the main 
theorem. The two lemmas above show that the sum on the left-hand side of 
l|A.14p equals 



J2 2f(^-^)detA„ = J2 



J2 det Xab , (A.16) 

i<J i<J \Ar\B=ir\J I 

which in turn equals the sum on the right-hand side of (|A.14p . 

detAAB. (A.17) 

Thus (|A.14p holds, and the theorem is proved. The final step from (jA.lGP to 
l|A.17P is justified by the observation that any given pair {A, B) of the type 
summed over in l|A.17P appears exactly once in the right-hand side of l|A.16p . 
namely for the sets / and J defined as follows. Let AI = An B, A' = A \ M, 
B' ^ B \ M, and let p > be the cardinality of the disjoint sets A' and B' 
(they are empty sets if p = 0) . Then define /' and J' by enumerating the 2p 
elements of A' U B' in the strictly interlacing order /' < J', and let I = M U I' 
and J = M U J' . Conversely, no other terms than these appear in the right 
hand side of (|A.16p . and it is therefore indeed equal to l|A.17p . 

Proof of Lemma \AA3[ The sets I < J and /' < J' (as in Definition I A. lip . with 

|/H|J|-fc, \I'\ = \J'\=p{I,J)^p, 

will be fixed throughout the proof, and for convenience we also introduce M = 
I n J and U = / U J, with \M\ = k — p and \U\ = k + p. We can assume that 
p > 0, since the case p = is trivial; it occurs when I = J, and then both sides 
of (|A.15P simply equal det A//. 

The set U \ M consists of the 2p numbers which belong alternatingly to /' 
and to J'. The sum (|A.15P runs over all pairs of sets {A, B) obtained by splitting 
these 2p numbers into two disjoint p-sets A' and B' in an arbitrary way and 
letting A = M U A' and B = M U B' . Write Q for this set; that is, Q denotes 
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the set of pairs (A, B) e C^"') x (I^;"') such that Au B = U and AnB = M. 
After expanding det Xab , we can then write the left-hand side of ijA.lSP as 

(-l)'^-'^aifc„(i)-'^a26„(2) ■ • -^Qfcb^Cfc)! (A. 18) 

where Sk is the group of permutations of {1,2,..., k}, and (—1)'^ denotes the 
sign of the permutation a. 

For each {{A, B), a) e QxSk, we let A' = A\M and B' = B\M, and set up a 
(cr-dependent) bijection between A' and B' as follows: a' e A' is paired up with 
b' € B' if and only if the product Xaib„^i-,Xa2b^^2) ■ ■ ■ ^a^b^^^) contains either the 
factor Xa'b' or a sequence of factors Xa'r, Xrs, • ■ • , Xtb' where r, s, . . . , t G M. 
Let us say that a' and b' are linked if they are paired up in this manner. A 
Hnked pair {a',b') e A' x B' will be called hostile if {a' ,b') belongs to /' x /' 
or J' X J', and friendly if (a', b') belongs to /' x J' or J' x I' . To each term in 
the sum (jA.lsp there will thus correspond p such Hnked pairs, and what we will 
show is that the terms containing at least one hostile pair will cancel out, and 
that the remaining terms (with all friendly pairs) will add up to the right-hand 
side of (|A15l) . 

Next we define what we mean by flipping a linked pair (a', b'). This means 
that, in the product Xa^b,(r,^a2b,(2) ■ ■ ■ ^a^b^^k) ■• those factors Xa'rXrs ■ ■ ■ Xty 
that Hnk a' to b' are replaced by Xyt ■ ■ ■ XsrXra', with all the indices in reversed 
order. (When the linking involves just a single factor Xa'b', fiipping means 
replacing it by Xb'a'-) Since the matrix X is symmetric, this does not change 
the value of the product, but it changes the way it is indexed. The number a' 
which used to be in the first slot (in Xa'r) is now in the second slot (in Xra'), 
and vice versa for b' . The connecting indices r, s, . . . ,t £ M do not contribute 
to any change in the indexing sets, since, for example, the r in Xa'r is moved 
from the second slot to the first, while the other r in Xrs is moved from the first 
to the second. The new product (the result of the fiipping) is therefore indexed 
by the sets 

{ai < ■■■ < flfc} 
{61 < • • • < h} 

respectively, and after reordering the factors so that the first indices come in 
ascending order, it can be written 

X~ T X~ T ... X~ r 

1105(1) £1205(2) afc05(fc) 

for some uniquely deterrnined permutation a G Sk- Flipping a given pair thus 
takes {{A,B),a) to {{A,B),(t). This operation is invertible, with inverse given 
by simply fiipping the same pair again, now viewed as a pair (6', a') e {{A)', {B)') 
Hnked via the indices t, . . . , s, r. Because of the symmetry of the matrix X, the 
term in (|A.18P corresponding to ((A,B),a) is equal to the term corresponding 



and 



A\{a'})u{b'}^:A^ 
(B\{b'})u{a'} 



B 
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to {{A, B),a), except possibly for a difference in sign, depending on whether the 
signs of a and a come out equal or not: 

We will show below that the permutation a has the same sign as a when a 
friendly pair is flipped, and the opposite sign when a hostile pair is flipped. 
Taking this for granted for the moment, divide the set 2x5^ into the two sets 
(Q X 5fc)hostiic, consisting of those {{A, B), a) for which at least one linked pair 
is hostile, and (Q x iSfe)friondiy, consisting of those {{A,B),a) for which all p 
linked pairs are friendly. The mapping "flip that out of all hostile pairs {a' ,b') 
for which min(a', b') is smallest" is an involution on (Q x 5fe) hostile that pairs up 
each term with a partner term that is equal except for having the opposite sign 
(since it is a hostile pair that is flipped). Consequently these terms cancel out, 
and the contribution from (Q x 5fc)hostiio to (|A.18P is zero. The sum therefore 
reduces to 

{-'^Y Xa^b,^r)Xa2b^^2) ■ ■ -Xakb^^k)- (A. 19) 

((yl,B),0-)e(Qx5fc)fric„dly 

Now equip the set (Q x 5fe)tricndiy with an equivalence relation; {{A,B),a) and 
{{A, B),a) are equivalent if one can go from one to another by flipping friendly 
pairs. Each equivalence class contains 2^ elements, since each of the p friendly 
pairs can belong to either /' x J' or J' x I' . Moreover, the terms corresponding 
to the elements in one equivalence class are all equal (including the sign, since 
only friendly pairs are flipped) , and each class has a "canonical" representative 
with all linked pairs belonging to /' x J', 

where the permutation a is uniquely determined by the equivalence class (and 
vice versa). Thus l|A.19p becomes 

2P J2 (-l)"^n,.<„^.2,.(2) • ■ -^.u.,., = 2PdetXjj, (A.20) 

which is what we wanted to prove. 

To flnish the proof, it now remains to demonstrate the rule that a has the 
same (opposite) sign as a when a friendly (hostile) pair is flipped. To this 
end, we will represent {{A,B),a) with a bipartite graph, with the numbers in 
U = ALi B (in increasing order) as nodes both on the left and on the right, 
and the left nodes ai G A connected by edges to the corresponding right nodes 
G B. The sign of a will then be equal to (—1)*^, where c is the crossing 
number of the graph. As an aid in explaining the ideas we will use the following 
example with U = [1,8], where the nodes in M = A C] B are marked with 
diamonds, and the nodes in A' and B' are marked with circles: 
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A = {2,3,4,5,6,8} 
= {2, 4, 5, 8} U {3, 6} 
= M U A' 



B = {1, 2,4,5,7,8} 
= {2,4,5,8}U{1,7} 
= MUB' 



Clearly, A' (J B' ^ {3, 6} U {1, 7} = {1, 3, 6, 7} ^ {i[ < j[ < i'^ < jj}, so that 
/' = {i'l, i'2} = {1, 6} and J' = ^2} = {3, 7}. Consequently, / = Af U /' = 
{1,2,4,5,6,8} and J = Af U J' = {2,3,4,5,6,7}. The chosen permutation 
is f7(123456) = 632415, where the notation means that a{l) — 6, cr(2) = 3, 
etc.; for example, the latter equality comes from the second smallest number 
02 in A being connected to the third smallest number 63 in B. There are 
9 crossings, so a is an odd permutation, and this graph therefore represents the 
term ~X2s.X^iXi2X5^XQiX%T , appearing with a minus sign in the sum IjA.lsp . 
The linked pairs (a', 5') £ A' x B' are (6,1) (directly linked) and (3,7) (linked 
via 4, 2, 8 e M). Both pairs are hostile, since (6, 1) S /' x /' and (3, 7) £ J' x J'. 

We will illustrate in detail what happens when the pair (3, 7) is flipped. 
The flip is effected by replacing the factors X34X42X28^87 by ^78-'^82-'^24-'^43 
and sorting the resulting product so that the first Jndices come in ascend- 
ing order; this gives X24X43^55^6i-'^78-'^82- Thus A = {2,4,5,6,7,8}, B = 
{1, 2, 3, 4, 5, 8}, and ^(123456) = 435162 (an even permutation). In terms of the 
graph, the nodes that are involved in the flip are, on both sides, {2, 3, 4, 7, 8} 
(the two nodes in the pair being flipped, plus the nodes linking them), and the 
edges involved are {34,42,28,87}, which get changed into {43,24,82,78}. In 
other words, the flip corresponds to this active subgraph being mirror reflected 
across the central vertical Hne. To understand how the process of reflection 
affects the crossing number, it can be broken down into two steps, as follows. 

On the left, node 7 is unoccupied to begin with, so we can change the edge 
87 to 77. This frees node 8 on the left, so that we can change the edge 28 to 88, 
which frees node 2 on the left. (Think of this edge as a rubber band connected 
at one end to node 8 on the right; we're disconnecting its other end from node 



32 



2 on the left and sliding it past all the other nodes down to node 8 on the left. 
Obviously the crossing number increases or decreases by one every time we slide 
past a node that has an edge attached to it.) Continuing Hke this, we get the 
result illustrated in Step 1 below; the edges changed are 87 — > 77, 28 — > 88, 
42 22, 34 44. 




^ ^ 0,,'^^ 

<e>- <> ^ 

Intermediate stage (after Step 1) Result of the flip (after Step 2) 

In Step 2, we work similarly on the right-hand side: node 3 is unoccupied to 
begin with, so we can change edge 44 to 43, and so on. The list of edge moves 
is 44 43, 22 24, 88 82, 77 78. (In the graph on the right we see 
that the crossing number after the flip is 8, verifying the claim that a is an even 
permutation.) 

We need to keep track of the changes in the crossing number caused by 
sliding active edges past nodes that have edges attached to them. This is most 
easily done by following the dotted lines in the figures, and counting whether 
the nodes that are marked (with circles and diamonds) are passed an even or an 
odd number of times. However, since the active subgraph simply gets reflected, 
the crossings among its edges will be the same before and after the flip, so we 
need in fact only count how many times we pass a passive marked node. (The 
passive nodes in the example are {1,5,6}.) 

If a passive node belonging to M is passed in Step 1, then it is passed the 
same number of times in Step 2 as well, since the nodes in M are marked both 
on the left and on the right. Therefore they do not affect the parity of the 
crossing number cither, and we can ignore the nodes marked with diamonds, 
and only look at the passive circled nodes (all the nodes in A' and B' except 
for the two active nodes that are being flipped). 

Passive nodes belonging to A' are counted only in Step 1 and passive nodes 
in B' only in Step 2; they get counted an odd number of times if they lie between 
the two flipped nodes (like node 6 in the example, counted once), and an even 
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number of times otherwise (like node 1, never counted). Consequently, what 
determines whether the parity of the crossing number changes is the number 
of nodes between the flipped ones that belong to A' U B' = /' U J'. And for a 
friendly pair, this number is even, while for a hostile pair, it is odd. 

This shows that the crossing number keeps its parity (so that (—1)'' = (— l)'') 
when a friendly pair is flipped, and the opposite when a hostile pair is flipped. 
The proof is flnally complete. □ 

B Verification of the Lax pair for peakons 

The purpose of this appendix is to carefully verify that the Lax pair formulation 
l|4.ip - l|4.2p of the Novikov equation really is valid for the class of distributional 
solutions that we are considering. This is not at all obvious, as should be clear 
from the computations below. 

B.l Preliminaries 

We will need to be more precise regarding the notation here than in the main 
text. A word of warning right away: our notation for derivatives here will differ 
from that used in the rest of the paper (where subscripts should be interpreted 
as distributional derivatives). 

To begin with, given n smooth functions x — Xk{t) such that xi{t) < ■ ■ ■ < 
Xn{t), let xo{t) = —CO and Xn+i{t) = +00, and let flk (for k — 0,. . . ,n) denote 
the region Xkit) < x < Xk+i{t) in the {x,t) plane. 

Our computations will deal with a class that we denote PC°°, consisting of 
piecewise smooth functions f{x,t) such that the restriction of / to each region 
Vlk is (the restriction to Slfc of) a smooth function f^^''{x,t) defined on an open 
neighbourhood of (so that /^'^-' and its partial derivatives make sense on 
the curves x — Xk{t)). For each fixed t, the function f{-,t) defines a regular 
distribution T/ in the class I?'(R), depending parametrically on t (and written 
Tf{t) where needed). After having made clear exactly what is meant, we will 
mostly be less strict, and write / instead of T/ for simpHcity. 

The values of / on the curves x = Xk{t) need not be defined; the function de- 
fines the same distribution T/ no matter what values are assigned to f{xk{t), t). 
But our assumptions imply that the left and right limits of / exist, and (sup- 
pressing the time dependence) they will be denoted by f{x^) := f'^^~^\xk) and 
/(x^) :— f'^^^Xk), respectively. The jump and the average of / at Xk will be 
denoted by 



[/(^fc)] f{4) - /(^n and {f{xk)) := ^ ^^'^ > ^ ^'^'^ > , (B.l) 



respectively. 



They satisfy the product rules 



(B.2) 
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We will use subscripts to denote partial derivatives in the classical sense, so 
that (for example) fx denotes the piecewise smooth function whose restriction 
to Slfe is given by df^'^^dx (and whose values at x — Xk{t) are in general 
undefined). On the other hand, will denote the distributional derivative, 
which in addition picks up Dirac delta contributions from jump discontinuities 
of / at the curves x = Xk{t). That is, D^Tf = Tf^ + J2k=i[fi^k)]Sxk, or, in 
less strict notation, 

n 

Dxf = f. + Y.[f(^k)]Sx,. (B.3) 

fc=i 

The time derivative Dt is defined as a limit in I?'(R), 

Tf(t + h)-Tf(t) , , 

D^Tfit) = lun^ 1 (B.4) 

and it commutes with by the continuity of Dx on V (R) . For our class PC°° 
of piecewise smooth functions, we have DtTf = Tf^ — X]fe=i ^fe [/(^fc)] '^a^t ) 
simply 

n 

fe=i 

where Xk = dxk/dt. We also note that f [x]^ (t) , t) — fx{xf {t),t) Xk{t) + 
ft{x^{t),t), which gives 

i[l{^u)\ = [L{^,)\ik+[ft{xk)l 
i{fixk))^{Uxk))ik + {ftixk)). 

B.2 The problem of multiplication 

If the function / is continuous at a; = Xk, then the Dirac delta at Xk can be 
multiplied by the corresponding distribution Tf according to the well-known 
formula 

TfSx^ fixk)6x,. (B.7) 

But below we will have to consider this product for functions in the class PC°° 
described above, where the value f{xk) is not defined. It will turn out that in 
the present context, the right thing to do is to use the average value of / at 
the jump, and thus define Tf 6x '■— {f{xk))Sx^. However, since we want this 
to be a consequence of the analysis, rather than an a priori assumption, we 
will, to begin with, just assign a hypothetical value f{xk) and use that value in 
l|B.7p . This assignment is justified in the present context, as we will see below. 
However, we are not claiming that this addresses any of the deeper issues; for 
example, this assignment does not respect the product structure of piecewise 
continuous functions. See [32l Ch. 5] for more information about the structural 
problems associated with any attempt to define a product of distributions in 
V'(R). 



35 



B.3 Distributional Lax pair 

Peakon solutions 



i(a:,t) = ^mfe(t)e-l---'=(*)l (B. 



belong to the piecewise smooth class PC°° . They are continuous and satisfy 



D^u = = ^mk sgn{xk ~ x)e 1^ ""=1, 
fc=i 

n n 

Dlu = D.^{ux) = Uxx + 'Y^[ux{xk)\ 5xk^u + ^(-2mfc) 5^^^ 

k=l k=l 

which implies 

n 

m:=u- Dlu = 2'^mk5xt,. (B.9) 

fc=i 

The Lax pair (|4.ip - (|4.2p will involve the functions u and D^u, as well as the 
purely singular distribution m. We will take V'l, V'2, '03 to be functions in PC^ , 
and separate the regular (function) part from the singular (Dirac delta) part. 
The formulation obtained in this way reads 

D^r* = A^' = A^, (B.IO) 

where * = (V"!, -02, "^s)*, 

/ n \ /O 1\ /O 1 0\ 

L = L + 2z^mfe4, MV, L= , iV= 1 ,(B.ll) 

\fc=i / \1 0/ \0 0/ 

and 

A = A-2ziY,'^ku{xkf5xAN, A={u/z -1/z^ -m^^/z . (B.12) 

\fc=l / \ —U^ u/z UUx J 

Note that (|B.10|) involves multiplying N'^ = (?/;2, V'3, 0) by Sx^, and some value 
ip2ixk) must be assigned in order for this to be well-defined (we will soon see that 
■03 must be continuous and therefore it is only ip2 that presents any problems). 

Theorem B.l. Provided that the product mip2 is defined using the average 
value ■ip2{xk) ■= {i'2{xk)) at the jumps, 

n 

mtp2 ■^2^mk{-ip2{xk))Sx^, (B.13) 

the following statement holds. With u and m given by (|B.8p - (|B.9p . and with ^' G 
PC°° , the Lax pair (|B.10p - i|B.12p satisfies the compatibility condition DfDx'^ = 
DxDt'^ if and only if the peakon ODEs (|3.4p are satisfied: Xk — u{xk)'^ and 
rhk = -nik u{xk){ux{xk)) ■ 
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Proof. For simplicity, we will write just ^ instead of X]fc=i- Identifying co- 
efficients of Sxf, in the two Lax equations (|B.10|) immediately gives [^'(x/c)] = 
2z mkN'^{xk) and — ifc [^"(0;/!;)] = —2zniku{xk)'^N'ii{xk), respectively. Thus, 
[V'3(a^fe)] — (in other words, is continuous) and Xk — u{xkY . Next we 
compute the derivatives of l|B.10[) : 

The regular part of (|B.10|) gives "^x = i^*, so that [A'^)x = ^i^* + AL'^ , 
and it is easily verified that LA — A^ + AL holds identically (since u^x = u). 
This implies that the regular parts of the two expressions above are equal, 
and the terms involving S'^^ are also equal since Xk — u{xk)'^. Therefore the 
compatibility condition Dt{Dx^) = Dx{Dt'^) reduces to an equality between 
the coefficients of Sxf, , 

-2zmku{xkfLN'^{xk) + 2zNf^{mk'i'{xk)) - [A^{xk)]. (B.14) 

Using the product rule l|B.2p . the expression for ['l'(xfe)] above, and [u2;(a;fe)] = 
— 2mfe, we find that the right-hand side of (|B.14p equals 

{A{xk))2zmkN^{xk)+[A{xk)]{'^{xk)) = 

2zmk I u/z I "^{xk) + 2mk I i/z 

' (*(.Tfc)). (B.15) 



u z 



The (3,2) entry — in the matrix in the first term will cancel against the whole 
first term on the left-hand side of l|B.14p , since the only nonzero entry of LN is 
(LiV)32 = 1. Thus (|B.14P is equivalent to 

TTik N-^{xk) + ruk Nj-^-^{xk) = 

(0 -u(^u^^ {u^^/z\ / u/z -1/2^ -2^ti^^/z\ 

u/z -1/z^ ^'(xfc) -f TOfe 1/2^ (*(2;fe))- (B.16) 
u/z \ -u/z 

To make it clear how the assumption (|B.13P enters the proof, we want to avoid 
assigning a value to i>2{xk) for as long as possible. Therefore we can't compute 
■^^{xk) quite yet. But (^{xk)) is well-defined, and its time derivative can be 
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computed using 'i'x = and = in (|B.6P : 

iV|(*(xfc)) - iV(L*(a;fc)) + Af(A*(xfc)) 

= iv(Lu(xfe)' + {A{xk))) {^{xk)) + N\ [Aixk)] [-^(xk)] 



(ujz -1/2^ -(u^)/z \ 
u/z u{u.) {n-k)) + \N[A{x,)]N2zm,n^,). 
00 ' 

A bit of manipulation using this result, as well as {i^3){xk) = i/^sixk), shows 
that the compatibility condition l|B.16P can be written as 

TOfciV^(^'(xfe) - (*(xfe))) + (riifc + mfcu(a;fc)(M^(xfe)))7V*(x) 

^mk(o u°/zo) h{xk)- ("iixk))) (B.17) 

The third row is zero, and the first two rows say that 

(mk + mku{xk){u.^{xk))^Tl'2{xk) = -"^fc^ (^/'2(a;fe) - {ip2{xk))^, 
(rrik + mku{xk){ux{xk))^'ip2,{xk) = \rnku{xk){i}2{xk) - (^'2(2;^))^. 

At this point we choose to assign i/'2(a;fe) := {ip2){xk), and then it is clear that 
l|B.17P is satisfied if and only if 

mfc = -mku{xk){ux{xk)) ■ 

□ 
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